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Abstract 

We investigate the vacuum in nonisentropic gas dynamics in one 
space variable, with the most general equation of states allowed by 
thermodynamics. We recall physical constraints on the equations of 
state and give explicit and easily checkable conditions under which vac- 
uums occur in the solution of the Riemann problem. We then present 
a class of models for which the Riemann problem admits unique global 
solutions without vacuums. 
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1 Introduction 

We consider the Euler equations of fluid dynamics in one space dimension, 

Pt + {pu)x = 0, 
{pu)t + ipu^ +p)x = 0, (1.1) 
(ipw^ + pE)t + (uiyu^ + pE+p))x = 0, 

describing conservation of mass, momentum and energy, respectively. Here p 
is the density, p is the pressure, u is the velocity, and E is the specific internal 
energy of the fluid: the equations are closed by specifying an equation of 
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state or constitutive relation, describing how the thermodynamic variables 
are related. This equation of state depends on the molecular structure of 
the fluid and is subject to physical constraints such as the Second Law of 
Thermodynamics. 

The Cauchy problem for (jl.ip is well understood when the initial data 
have small total variation [Hj \T\ \16\\7\. but little is known for solutions with 
large data \18\ I19j. One of the main difficulties is the possible occurrence of 
vacuum. Several authors have studied the vacuum in isentropic gas dynam- 
ics (obtained by dropping the energy equation) [101 [El ESI [20] , and more 
recently nonisentropic polytropic ideal (7-law) gases [3llll[5]. 

Our goal here is to describe in detail which gases admit vacuums in the 
solution. It is well known that 7-law gases require vacuum in order to solve 
the global Riemann problem [16], whereas an isothermal gas (7 = 1) does 
not. Also, a better existence theory is available for the isothermal system, 
due to a degenerate wave curve structure [14J, but there is no natural 3x3 
physically consistent analogue of the isothermal system. In ^17j, Temple 
considers a class of 3 x 3 constitutive relations with simplified structure, 
but this system violates some thermodynamic constraints, given in [11^ I13j. 
Our intention is to present a class of physically consistent constitutive laws 
which do not admit vacuum, thus removing the issue of vacuum for these 
equations of state, in order to focus more fully on the effects of nonlinear 
wave interactions. 

We begin by collecting all the physical conditions that restrict the equa- 
tion of state. We recall the solution of the Riemann problem and give a 
necessary and sufhcient condition for the occurrence of vacuums in the gen- 
eral solution of the Riemann problem. By exhibiting specific examples, we 
describe a class of constitutive laws which satisfy all our physical constraints 
and for which the Riemann problem does not contain a vacuum state. We 
then present the simplest class of such equations of state, namely 



where 2 < c ^ 1- Although it is implied by quoted results, we explicitly 
solve the global Riemann problem for this pressure law, without use of the 
vacuum state. It is expected that this will inform such a study of wave 
interactions, which is the subject of the authors' ongoing research. 
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2 Thermodynamic constraints 

To set notation, we describe the thermodynamic constraints. The ther- 
modynamic properties of a fluid are embodied in the constitutive relation 
E = E{t,S), where r = 1/p is the specific volume and S is the specific 
entropy. See [13] for a physical discussion of these constraints, and |15j for 
a detailed mathematical analysis. 

First, we make the smoothness assumption, 

E = E{t, S) G C^, for r G M+ and 5 G M, (2.1) 

which is true for most fluids. 

We require the fluid to satisfy the Second Law of Thermodynamics, 
which asserts that 

clE = TdS -pdr, (2.2) 
where T is the temperature. This in turn implies that 

Es = T, and Er = -p. (2.3) 

We assume the standard thermodynamic constraints: specific vol- 
ume T, pressure p and temperature T satisfy 

r > 0, p>0 and T ^ 0; 

so that by (fOjl . 

Er <0 and Es ^ 0. 
We assume "stability of matter" , which asserts that the energy is finite, 

£^00 = lim E < 00, 

and without loss of generality we take 

Eoo = 0, so that E{t, S) > (2.4) 

for all r > and S, see JTT\. 

Next, we assume the thermodynamic stability constraint that the 
energy E = E{t, S) be jointly convex. 

Err = -Pr > 0, Egs = Tg ^ 0, (2.5) 

while also 

Err- Ess 1^ Els and ErS = -ps^Q- (2.6) 
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According to [13j, thermodynamic stability yields (|2.5p and ()2.6p . but our 
discussion requires only (j2.5p . Equation (j2.5p in turn implies that the sys- 
tem is strictly hyperbolic away from vacuum. The condition ps > states 
that the material expands upon heating at constant pressure. We assume 
nonstrict inequality for ps to include isentropic gas dynamics, for which 
PS = 0. 

Our final condition is an energy condition, which states that if the 
pressure p{t, S) is specified, then the energy E is well defined: that is. 



/oo 
pir',S) dr', 



where we have used (|2.3p and ()2.4p . That is, we require that, for all S, 

\{T,S)dT= [' dp < +00, (2.7) 
Jo P 

where p is defined by 

p{p,S)^p{l/p,S)=p{T,S). (2.8) 

The energy condition (12. 7p imposes growth conditions on p{t,S), or 
equivalently restricts the pressure p near vacuum, namely ^3(0+, S) = 0, 
and by rHospital's rule, 

lim pp{p, S) = lim = 0. (2.9) 

p->-0 p-i-O p 

Note that our conditions alone are not sufficient to conclude uniqueness 
of solutions to the Riemann problem: uniqueness is assured if and only if 
Smith's medium condition, that is 

where p is regarded as p{t,E), is satisfied, see |15j . 



3 Vacuum in the solution of Riemann problems 

We wish to investigate the circumstances in which a vacuum appears in 
the solution to a Riemann problem. We briefly recall the solution of the 
Riemann problem; see |16j . 
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3.1 Riemann problem 

We begin by calculating the simple (rarefaction) wave curves. For smooth 
solutions, we replace the third (energy) equation of (jl.ip by the entropy 
equation 

St + uSx = 0, 

and use {p, u, S) as the state variables. It is routine to calculate the eigen- 
system after writing (jl.ip in quasilinear form. The eigenvalues of (jl.ip are 

Xi = u — c, X2 = u, A3 = + c, 

and these are the wavespeeds of the backward, middle and forward waves, 
respectively, and 

c = c{p,S) := (3.1) 

is the speed of sound. As is well known, the forward and backward waves 
are genuinely nonlinear and the middle waves linearly degenerate. The cor- 
responding eigenvectors are 





'^i = I -c I , r2 = I U I , rs 

-Pp 

It follows that the equation of a backward simple wave is 

u-ui = R{pi,S)-Rip,S), S = Si, (3.2) 
where the subscript / refers to the left state of the wave, and we define 

R(p,S)= [^^^^dr= fy^^dr, for p ^ 0. (3.3) 
Ji r Ji r 

The equation of a forward simple wave curve is 

Ur — U = R{pr, S) — R{p, S), S = Sr, (3-4) 

where the subscript r refers to the right state of the wave. 

Next we calculate the shock curves: these are described by the Rankine- 
Hugoniot conditions, 

= [pu], 

C[pu] = [pu^+p], (3.5) 
C[^pu^ + pE] = [ui^pu^ + pE + p)], 
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where ^ is the shock speed and the brackets denote the jump in a quantity 
across the shock. We simphfy these as follows: the first equation can be 
written 

- C - ^n. + = 0, (3.6) 
Pi Pi 

and, recalling that r = using this in the second equation and simplifying 
yields 

[p][T]=i[pu][T]-W][T] = -[u]\ 

Next, denoting the average of a quantity by 5 = , manipulating the 
first two equations of (j3.5p yields 

_ _ _ _ b] 

ip[u\= pu[u] + \p\, so ip-pu = Y^. (3.7) 

The third equation of (j3.5p gives, after simplifying, 

i'p[\u'^ + E]= lm[\u'^ + E] + [up]. 

Using (j3.7p and again simplifying, we finally obtain 

[E]+p[t] = Q, (3.8) 

which is the Hugoniot curve for shocks. We conclude that this describes 
the shock curve fully: first, solve (j3.8p to find the relation between p and 5, 
then use 

N = -VIpIF], (3.9) 

obtained from the entropy condition [9l[6], to resolve u, and finally use (j3.7p 
to determine ^. 

Recall that an entropy condition is required to choose admissible shocks 
and thus obtain uniqueness of Riemann solutions [Bl |9]. This condition 
states that pressure (and thus also density) is bigger behind the shock, and 
leads to the negative square root in ()3.9I) above. It follows similarly that the 
density behind a (forward or backward) rarefaction wave is smaller than the 
density ahead of the wave. Also, this implies that the state behind a shock 
cannot be the vacuum state. 

It is routine to describe contact discontinuities using (j3.5p : namely, sub- 
stitute [u] = in directly, to obtain 

[u] =0, [p] = and ^ = u = u. (3.10) 
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If a contact discontinuity is adjacent to a vacuum, then we combine 
the contact discontinuity and the vacuum region into a new vacuum region, 
cahed a non-isentropic vacuum on which the entropy density pS vanishes. 
The left and right hand hmits of S on the left and right boundaries of such 
a vacuum region are different. It follows that, if the vacuum is involved in 
the solution of the Riemann problem, it can only be generated between two 
outgoing rarefaction waves, see [121 [20] 

3.2 Vacuum condition 

Lemma 3.1. The vacuum state exists in the solution of Riemann problems 
of if and only if, for some S, 

R{0+,S)>-oo, (3.11) 

where R is defined in i3. 

Proof. We first prove that if the vacuum state exists in the solutions of 
Riemann problems then ()3.1ip is satisfied. Recall that vacuum state only 
appears between two rarefaction waves. If a Riemann solution consists of 
forward and backward rarefactions, it follows that the velocity u is monotone 
increasing as a function of x [6]. Parameterizing the forward wave by p, it 
follows from (|3.4p that 

Ur — U = R{pr, Sr) — R{p, Sr)', 

now since u ^ ui, we get the uniform bound 

R{p, Sr) ^ Ui — Ur + R{pr, Sr), 

and allowing p — > implies ()3.1ip . 

Now suppose (j3.1ip holds for some S. We claim that the Riemann 
problem with data Ui = (1,0,5"), and Ur = {l,Ur,S) has a vacuum in 
the solution whenever 

Ur > -2R{0+,S). 

To see this, assume Ur > and resolve the Riemann problem into backward 
and forward rarefactions using (j3.2p and p.4p . to get 

Um = R{l,S) - R{pm,S) = -R{pm,S) and 

Ur-Um = R{1, S) - R{pm, S) = -R{pm, S), 



7 



with no contact as Sr = Si = S. Adding, we must solve 



Ur = -2R{pm,S), 

and so if Ur > —2R{0+,S), no such pm can be found and a vacuum is 
required to solve the Riemann problem; see also \16\ I20j. □ 

We now introduce an easily checkable condition which implies (j3.1ip . 
so is a sufficient condition for existence of Riemann solutions with vacuum. 
This pressure near vacuum condition describes the rate at which p ^ 0: 
for some value of S, there exist positive numbers Eq, uq and Mq, such that, 

Pp{p,S) !^ Mop""" whenever /jG(0,eo). (3.12) 

Note that we require this condition at only one S: by continuity, we would 
generally expect the condition to hold in an open set of S values. Note also 
that polytropic ideal gases satisfy ()3.12p . ao being given by the adiabatic 
exponent 7 > 0. 



Theorem 3.2. The energy condition (2.1) and pressure near vacuum con- 
dition 1^3. together imply that vacuums exist in the solution of some Rie- 
mann problems. 

Proof. According to (j2.9p . the energy condition implies that — > as 
p — 7- 0, so ()3.12p makes sense. It then suffices by Lemma l3. II to show that, 
for S given by (j3.12p . equation (j3.1ip is satisfied. From (|3.3p . for p < eq, we 
have 



R{eo,S)-Rip,S)= Ty^^^dr 



^ / ^/M'Qr-^+^'■'/'^ dr 



ao 

and taking the limit p — >■ gives the required lower bound for i?(0+, S). □ 



4 Gas dynamics without vacuum 

We now write down a class of constitutive laws for gases which do not admit 
vacuums in the solution of the Riemann problem. These gases satisfy all the 
constraints of Section 2, but do not satisfy the vacuum condition (j3.1ip . 
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By (j3.12p . there is a vacuum for any 7-law gas (p ^ p"') with 7 > 1, but 
no vacuum for an isothermal gas (7 = 1). We thus look for presssure laws 
between these cases, which restricts our equation of state. 

It is convenient to work with separable energies, that is, energies of the 
form E{p,S) = f{S)g{p). Specifically, we consider energies given by the 
expression 

r+00 

Eir,S) = f{S) k\y)dy, (4.1) 

Jln(r+1) 

where / and k are functions satisfying certain conditions described below. 
We assume that f{S) is C^, positive, increasing and convex, 

f{S)>0, nS)>0, r(5)^0, (4.2) 

with limits 

lim f{S) = 00 and lim f{S) = 0. (4.3) 

We assume k{y) is a C^, convex, decreasing function defined on (0, 00), 

k{y)>0, k'{y)<0, k"{y)>0, (4.4) 
subject to the growth limits 

/+00 r+oo 
k'^{y) dy < +00 and / k{y) dy = +00, (4.5) 

describing growth near vacuum, r — )• 00, and limits 

k'{y) 

kiy) — 7> 00 and -r-r-\ — ^ ^.s y — )• 0, (4-6) 

Kv) 

describing infinite density. 

Theorem 4.1. The above conditions imply that for gases having energy 
the Riemann problem has a unique global solution which does not 
admit the vacuum state. Moreover, such gases satisfy all constraints of 
Section 2 with the exception of [27t 



Proof By ([H]) and ([23]), we have 

piT,S) = h\ln z)f{S), (4.7) 

where we have set z = r + 1. It is easy to check that the standard thermo- 
dynamic constraints, stability of matter, and thermodynamic stability (j2.5p 
are satisfied. The energy condition follows from the first equation in 
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We now check that the vacuum condition (|3.1ip fails for this equation 
of state. We calculate 

fc^(lnz)-2fc(lnz)fcXlnz) ,,,, 

Pr[T,b) = 2 < 0, 

and, since p = l/r, we write 

PpiPiS) = -Pr{T,S)T'^, 

and so, by ()3.3p . 

R{p,S) = -j V-Prir,S) dr 

, /-InCl+l/p) , 

= -\/7(50 / -2fcHA:'(w;) dw, (4.8) 

Jin 2 

where w = Inz = ln(l + r). It now follows from (14. 4p . (I4.5P that 

, /-InCl+l/p) 

i?(/9, 5) < -Vf{S) / A;H ^ -oo 

J In 2 

as p — )• for all S, so that (|3.1ip is never satisfied. 

Existence and uniqueness of Riemann solutions with arbitrary data now 
follows from Smith's medium condition (j2.10p . see [15j . Smith has an extra 
assumption, namely, he requires 

lim E(p, r) = 0, (4.9) 

T-S>0+ 

when E is described as -E = E(p,t). Because our energy is separable, we 
have 

Since z — >■ 1 as r ^ 0, for small r, we write 

By our assumptions on k, (|4.9p follows if we show that the limit 

k'^(y)dy 

, ^0 as e ^ 0. 
k^{e) 
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This in turn follows from (j4.6p and 1 'Hospital's rule. 

To check the medium condition (j2.10p . we use (j4.10p to write 

, ^ k^Clnz) 

p{t,E)=E ^ ^ 



where z = t + 1. The medium condition ()2.10p is 

2k{\nz)k'{\nz) kHy) dy - k\\nz) £^ k\y) dy + fc^(lnz) l_ 
(zjZk\y)dy? ^2E' 

which simphfies to 

N{z) := -k\z) + 2k' {z) / k\y) dy - k{z) / k\y) dy < 0, 



which must hold for all z > 1. 

By ()4.5p . k{z) ^ as z ^ oo, which in turn implies that N{z) ^ as 
z —7- OO. Furthermore, by 



-1 poo pc 

N'{z) = --k'e + 2k" k^-k' k^ + k^>0, 

2 Jz Jz 



so N is increasing with limit 0, and thus N[z) < for all z. □ 

We remark that our conditions do not suffice to prove the convexity of 
E{t,S), namely the first condition of (j2.6p . 



Err ■ Ess ^ Els- 

This condition is easily seen to be implied by the dual assumptions that / 
is log-convex, 

f{S)nS)^ f'\S), (4.11) 
and k satisfies the condition 

poo 

{k{z) - 2k'{z)) J k^{y) dy ^ k^{z). (4.12) 

for all z. 
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5 Concrete example 



We now present the simplest equation of state which does not allow for a 
vacuum in the solution of the Riemann problem. Ideal polytropic (7-law) 
gases with adiabatic constant 7 > 1 admit vacuums, while an isothermal gas 
(7 = 1) does not, but the isothermal gas does not satisfy the finite energy 
stability of matter condition. We thus consider equations of state that fall 
between these two cases. Since (|3.12|) implies existence of vacuums, this 
restricts the equation of state to those for which ()3.12p fails. 
This class of gases we present here are polytropic, satisfying 

E = CrT, 

but do not satisfy the ideal gas law pr = RT. For a polytropic gas of the 
form (|4.ip . (j2.3p implies that we must have 

f{S) = Koe^'^\ 

and this trivally satisfies conditions (|4.2p . (j4.3p and (j4.1ip . 

We choose k = k{z) as simple as possible so that properties (14. 4p and 
(j4.5p hold, namely 

k{z)=z~", for some a € {^,1]. (5.1) 

It is then easy to check that ([331), (jlSl), <SM and hold. 
With these choices, ()4.1|) becomes 

E = ^^-^{Hr + l))'-'% (5.2) 

and (|4.7p becomes 

p=^— -(ln(r + l))-2^ (5.3) 
r + 1 

It follows from Theorem 14.11 that the Riemann problem has a unique 
solution without vacuum. In fact, Smith's strong condition holds, namely 

This is easily seen by eliminating e^^'^^ from (|5.2p and (|5.3p . to get 

(r + 1) ln(T + 1) 
Eir,p)=p ^^-^ , 

and differentiating. 

Although we have an abstract proof of existence and uniqueness of Rie- 
mann solutions, we find it instructive to prove this directly. 
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Lemma 5.1. There is a unique solution to the Riemann problem for il.l]) 
with equation of state i5.^) . which does not include the vacuum, for arbitrary 
Riemann data. 

Proof. Using (j4.8p and simplifying, we have 

/■ln(r+l) , 

i?(r,5) = -v^e^/2c. / ^-'^Ji+2^dw. (5.4) 

J In 2 ^ 

Following [2], we now make convenient changes of variables. First, set 
4> := ln(r + 1) > 0, so that T + l = e^, 

and define 

h := e-"^/2c. ^(^^ ^ - /'^ w-" Jl + ^ dw. 

J In 2 ^ 

It is clear that p, r, (j) ^ £^re equivalent coordinates, with and r 
decreasing and /i increasing as functions of p, and since a G 1]; we have 

lim h = —oo and lim h = oo. (5-5) 

Next, define 

so that {(j), m) can be used in place of (r, 5) or (p, S). 
It follows from IK2\j and (EHl) that 



E = - — ^0m^ and p = i^o e rn^ , (5-6) 
2(j — 1 

while also 

gS/2c, ^ ^ g^j^^ R = hm(j)''. 
The simple wave curves p.2p . ()3.4p are described by 

Ur-Ul = [ha - hb)ma(pa, = T^' (^-'j 

rua (pl 

where the subscripts denote the behind, ahead, right and left states respec- 
tively, and (f) = (j){h). For rarefaction waves, the sound speed c decreases 
from front to back, so these are characterized by ha > h^. Similarly, by 
p.lOp . a contact discontinuity is described by 

ur = uu ^ = e(<^^-<^')/2^ (5.8) 
mi 
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It remains to calculate the shock curves. Using ()5.6p in ()3.8p . we get 

{(t>bmb - (pama'') + Kq (e'^i' - e"^") = 0, 



2a -I 
which becomes 



m, 



' ' 5- 4). (5.9) 



Here, 4>a > (pb, or equivalently ha < hb, since the sound speed c is greater 
behind the shock O |9]. It is clear that the function / in (j5.9p makes sense 
only if the function inside the square root is nonnegative. To check that this 
holds, consider the function 

for X > y > 0. Denote the numerator and denominator of q by 

qn{x,y) := x H e^~^ and 

2(7 — 1 2 2 

Qdix,y) ■■=^,y + \-r^-'^ 

respectively. It is immediate that for x > 0, 

qn{x,x) > 0, qn{x,0) > 0, qd{x,x) > 0, qd{x,0) < 0. 

Since qd{x, y) is strictly increasing with respect to y, for each x, qd{x^ y) = 
has a unique solution < y = ^{x) < x. On the other hand, qn{x, y) is 
strictly decreasing with respect to y, so qn{x, y) > for x > y > 0. Hence, 
for fixed x, 

q{x,y) > as long as ¥?(x) < y < x. 



while also 



lim (/(x, y) = +00. 



It follows that, for fixed (pa, < (p{(j)a) < (pa and the shock curve is 
parameterized by (pb G {^{(pa), (pa), with 

lim f{cPa,<Pb) = +oo. (5.10) 

<t>b^'fii<t>a) + 
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Moreover, f{(j)a,4'b) increases with respect to (j)a and decreases with respect 
to (pb. Since qd{x,y) = qn{y,x), we also have 



f{4>b,(t>a) 

Next, in these coordinates, (j3.9p is 

[u] = -\fK^\J {e-'^brrib'^ _ e-'^^ma'^){e't>^ - e't">) 

= -y/K^J {e^--<l>i>^^ - 1)(1 - e^b-<l>a) rua. (5.11) 

Defining 

ff(0a,4) := -v^^(e<^-<^<'/2(</.„,/.fe) - 1)(1 - e't'b-'t'^), (5.12) 
it is easy to check that 

9{4>a,4>b)ma = g{4>b,(l)a)mb, 

provided (j5.9p holds, and g{(j)a, 4>b) decreases with respect to and increases 
with respect to 4)b- By (|5.10p . since (pa > 4>b-, 

lim g{(t)a, h) = -oo. (5.13) 
To express the composite wave curves, we define 

G(/^„/.,):=^^'^-'^^)' ha<hb<h{^{4>a)), 

I {ha - h)(t)^, ha ^ hb, 



F{ha,hb) :-- 



f{4>a,(t>b), ha < hb < h{(p{<f>a)), 
ra/rb, ha ^ hb, 



with = 4>(K), and '^{x) is defined by qd{x,ip{x)) = 0. Both F{ha,hb) and 
G{ha,hb) are defined on the region {hb < ip{ha)}, with F increasing and G 
decreasing in hb for fixed ha- Moreover, by ()5.5p . we have the limits 

lim F{ha,hb)=0 and lim G{ha, hb) = +oo, (5-14) 

hf,— oo hj,— oo 

for any ha fixed, and by (|5.10p and (|5.13p . 

lim F{ha,hb) = +OG and lim G{ha, hb) = —oo. (5.15) 

ht^h(!fi((f>a)) hi,^h{tp((j)a)) 
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Combining ()5.7p and ()5.9p . ()5.1ip . we describe the composite 1- and 
3-wave curves by 

Ur - Ul = G{ha,hb)ma, — =F{ha,h), (5.16) 

while the 2-waves are given by (jS.Sp . namely 

rrir <t>r-4>i 

Ur = Ui, = 62 . 

mi 

We now solve the Riemann problem by resolving the intermediate states. 
We use subscripts L, 1, 2 and R to denote the left, intermediate, and right 
states, respectively. We set 

A:="^-"^ and B:=^, 
mi rriL 

and use (|5.16p and (|5.8p to describe the waves, and eliminate u, to get the 
equations 

G{hL,hi) + BG{hR,h2) = A, and (5.17) 
FQiL, hi) e-^^l'' = B FQiR, /is) e''^^/^^ (5.18) 

and we must solve for hi and /i2- 

For any fixed /iq, the function FiJiQ^h) e~'^^^^/'^ is increasing in /i, and 
by (15.14p . (I5.15p . it has range (0, oo). Thus there exists a unique (increasing 
and C^) function of /12, such that 

hi = G(/i2;-B,/iL,/ii?) 

if and only if ()5.18p holds. Substituting into (j5.17p gives the equation 

G{hLMB, hL, Hr, h2)) + B G{hR, h2) = A, 

and we must solve for /12. Since G{hQ,h) is decreasing with h, and by 
(j5.14p . (|5.15p has range (—00,00), there is a unique solution /i2 < ^{4>r) of 
this equation. Note that we also have hi < ip{(f)i). Finally, we use (15.16P to 
fully resolve the intermediate states, and the construction of the Riemann 
solution is complete. □ 



16 



References 



[1] A. Bressan, Hyperbolic Systems of Conservation laws: The One- 
dimmensional Cauchy Problem, Oxford Lecture Ser. Math. Appl. 20, 
Oxford Univ. Press, Oxford, 2000. 

[2] Ceng Chen, Formation of singularity and smooth wave propagation for 
the compressible Euler equations, to appear in J. Hyp. Diff. Eq.. 

[3] G. Chen, E. Endres and K. H. Jenssen, Pairwise wave interactions in 
ideal polytropic gases, to appear in Arch. Rat. Mech. Anal. 

[4] G. Chen and R. Young, Smooth solutions and singularity formation for 
the inhom,ogeneous nonlinear wave equation, J. Diff. Eq., 252(2012), 
2580-2595. 

[5] G. Chen and R. Young, Shock formation and exact solutions for the 
compressible Euler equations with piecewise smooth entropy, submitted. 

[6] R. Courant and K. O. Friedrichs, Supersonic Flow and Shock Waves, 
Wiley-Interscience, New York, 1948. 

[7] C. M. Dafermos, Hyperbolic Conservations laws in Continuum Physics, 
Springer- Verlag, Heidelberg, 2000. 

[8] J. Glimm, Solutions in the large for nonlinear hyperbolic systems of 
equations. Comm. Pure Appl. Math., 18(1965), 697-715. 

[9] P. D. Lax, Hyperbolic systems of conservation laws, U, Comm. Pure 
Appl. Math., 10(1957), 537-566. 

[10] T.-T. Li and Y. Zhao, Vacuum problems for the systems of one- 
dimensional isentropic flow (in Chinese), Chinese Quarterly Journal 
of Math., 1(1986), 41-46. 

[11] E. Lieb, The stability of matter. Rev. Mod. Phys. 48(1976), 553. 

[12] T.-P. Liu and J. Smoller, On the vacuum state for the isentropic gas 
dynamics equations. Adv. in Appl. Math., 1(1980), 345-359. 

[13] R. Menikoff, B. J. Plohr, The Riemann problem for fluid flow of real 
materials. Reviews of Modern Physics, 61:1(1989), 75-130. 

[14] T. Nishida, Global solutions for an initial boundary value problem of a 
quasilinear hyperbolic system, Proc. Japan Acad., 44(1968), 642-646. 



17 



[15] R. Smith, The Riemann problem in gas dynamics, Trans. Am. Math. 
Soc, 249(1979), 1-50. 

[16] J. SmoUer, Shock waves and reaction- diffusion equations, Springer- 
Verlag, New York, 1982. 

[17] B. Temple, Solutions in the large for the nonlinear hyperbolic conser- 
vation laws of gas dynamics, J. Diff. Eqns., 41(1981), 96-161. 

[18] B. Temple and R. Young, The large time stability of sound waves. 
Comm. Math. Phys., 179(1996), 417-466. 

[19] Blake Temple and Robin Young, A paradigm for time-periodic sound 
wave propagation in the compressible Euler equations. Methods and 
Appls of Analysis, 16(2009), 341-364. 

[20] R. Young, The p-system IT. The vacuum, Evolution Equations Banach 
Center Publications 60(2003), 237-252. 



18 



